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Abstract. Stationary waves in a superfluid magnetoexciton gas in = 1 quantum 
Hall bilayers are considered. The waves are induced by counter-propagating electrical 
currents that flow in a system with a point obstacle. It is shown that stationary 
waves can emerge only in imbalanced bilayers in a certain diapason of currents. It 
is found that the stationary wave pattern is modified qualitatively under a variation 
of the ratio of the interlayer distance to the magnetic length djl. The advantages of 
use graphene-dielectric-graphene sandwiches for the observation of stationary waves 
are discussed. We determine the range of parameters (the dielectric constant of the 
layer that separates two graphene layers and the ratio dj t) for which the state with 
superfluid magnetoexcitons can be realized in such sandwiches. Typical stationary 
wave patterns are presented as density plots. 
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1. Introduction 

Realization of superfluidity of excitons in semiconductor heterostructures is considered 
as a challenge problem. Exciton mass is smaller or of order of the electron mass and the 
temperature of the transition into a superfluid state can be quite large. Considerable 
attention has been given to the superfluidity of spatially indirect excitons in bilayer 
electron structures. A flow of such excitons results in an appearance of counterflow 
electrical currents in the layers. They can provide feeding of the load located at one end 
of the system from the source located at the other end, and superfluidity of spatially 
indirect excitons can be considered as a kind of unconventional superconductivity. The 
idea goes back to the papers lU |2]. This idea was reincarnated in the quantum Hall 
physics of bilayer electron systems in Refs. [31 HI [5] . In a quantum Hall bilayer with 
the total fllling factor of Landau levels i/^ = 1 the state with a spontaneous interlayer 
phase coherence may emerge (see P [71 El El [10] as a review). The many-particle wave 
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function for such a state describes the pairing of electrons from one layer and holes from 
the other layer and its structure is analogues to the BCS wave function. Experimental 
studies of transport properties of quantum Hall bilayers [HI [121 HSl [lH [IS] confirm in 
some part the theoretical prediction. The systems demonstrate vanishing of the Hall 
resistance and exponential increase of the counterflow conductivity under lowering of 
temperature. But the conductivity remains finite that can be accounted for unbound 
vortices [l6l|l7l[T8]. 

Another hallmark of superfluidity is that an obstacle does not disturb a flowing 
superfluid at small flow velocity. If the flow velocity exceeds the critical one, the 
obstacle produces excitations that are in rest relative to it. They are the stationary 
waves - the waves with the time independent phase in each space point (in the reference 
frame connected with the obstacle). Stationary waves may emerge in a medium which 
excitation spectrum has a dispersion. The well-known example is ship waves - the 
waves produces by a ship on the water surface. The spectrum of deep water waves has 
negative dispersion and ship waves on a deep water are located behind the floating ship. 
The Bogolyubov spectrum has positive dispersion and stationary waves in a superfluid 
weakly non-ideal Bose gas are located outside the Mach cone: the waves appear ahead 
the obstacle, and do not appear strictly behind it [20j. 

Stationary waves have been recently observed in an exciton-polariton gas [21] . Since 
the exciton-polariton mass is much smaller than the electron mass, a two-dimensional 
exciton-polariton gas in microcavities is considered as a candidate for a high-temperature 
superfluid [22l |23l [2lj. The observation [21j consists in that at rather high exciton- 
polariton concentration and small flow velocity the waves are not excited. But if the 
flow velocity exceeds the speed of sound stationary waves emerge outside the Mach cone. 

In this paper we consider the possibility of observation of stationary waves in 
quantum Hall bilayers in a state with a spontaneous interlayer phase coherence. In 
Sec. [2] we use the kinematic approach [19] based on the analysis of the excitation 
spectrum and determine the range of the fllling factor imbalance and the counterflow 
currents where stationary waves may emerge. It is found that the stationary wave 
pattern is modifled qualitatively under variation of the ratio of the interlayer distance 
d to the magnetic length i. In Sec. [3] we present the description of stationary waves in 
ut = 1 quantum Hall bilayers in the dynamical approach [19]. The approach allows to 
compute the amphtudes and phases of stationary wave in each space point. The results 
are presented as density plots. In Sec. |4]we discuss advantages of use graphene bilayers 
for the observation of stationary waves. 

2. The excitation spectrum and conditions for the emergence of stationary 
waves in quantum Hall bilayers 

In this section we consider under what conditions stationary waves can be excited in a 
superfluid magnetoexciton gas in bilayers. 

We use the kinematic approach [19] that is based on the analysis of the excitation 
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spectrum. The spectrum at general direction of the wave vector relative to the direction 
of counterflow currents and general filling factor imbalance was obtained in [251 126] 
basing on the approach [27]. Here we reproduce in short the derivation p6]. 

The state with the spontaneous interlayer phase coherence is described by the 
many-particle wave function [6] 

1^) = n + T<^-««^0 

X ^ ^ 

Here we use the form of ([1]) given in Ref. [26]. It corresponds to a coordinate 
system with the x axis directed at some angle to the counterflow current direction. 
The quantities a\-^^ ai^x are the creation and annihilation operators for the electrons 
in zeroth Landau level (the lowest Landau level approximation is used), % is the 
layer index, and X is the guiding center coordinate. The vector-potential is chosen 
in the form A = (0, -Bx, 0) that corresponds to the single particle wave function 
'?/'x(r) = (l/vr-'^/'^£)e*'^^/^%xp[— (x — X)^/2£^]. The variable dx can be presented as 
a sum of the mean value and the fluctuating part: = + ^x- The mean value Qq 
is connected with the filling factors by the relation cos^o = — ^2, where z/j = li^nii^ 
is the filling factor for the i-th layer, and nj is the electron concentration. The phase 
contains the regular Q^X and the fluctuating (px parts. 

The counterflow supercurrents are directed along the vector Q = {Qx.Qy)- Their 
values are given by the equation 

ScF = (2) 

The quantity jcF is defined as the density of the supercurrent in the layer 1 (ji = 
~j2 = jcF), Eq = (\Efo|-f^c|^o) is the mean-field energy for the Hamiltonian of Coulomb 
interaction in zeroth Landau level in the state ([I]) at 6*^ = <px = 0, and S is the area of 
the system. 

The explicit expression for Eq reads as 

^0 = ^ [cos^ eo[Ho - FsiO)] - sin^ eoFn{Q)] , (3) 
where the quantities 

Es{D){q) = ^ pMP(l^'^)Vs(D){p)e~'—dp (4) 

{Jo{x) is the Bessel function) describe the contribution of the intralayer (S) and the 
interlayer (D) exchange Coulomb interaction, and the quantity 

the contribution of the direct Coulomb interaction. Vs{q) = 2TTe'^/eq and Voiq) = 
2'Ke^e~'^'^ / eq are the Fourier components of the interaction, and e is the dielectric 
constant. 
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Q 

-{im,{-q)K,^{q)(p{q) + c.c.)], (6) 
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where 



1 /27r/2 ^ 



2V S ^ 



The components of the matrix K are defined as 



(l=qix 



where 



K,,(q) = i/(q, Q) - F^dql) + F^(|Q|) 

■Fn{\q-< 

i^D(|q+Q|) + FD(|q-Q|) 



+ (^.(|Q|)-M^L±M±MM)eot^,, 



-ft'^^(q) = sin^ 6', 



i^D(IQI) 
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^^v'lq) = -cos^o ^ , (8) 



and 



^(q, Q) = ^ - (lq >< Ql^')] e""^- (9) 

The quantities hmz{q) and ^{—q) are the conjugated variables. It can be checked 
by computing the commutator of the operators that correspond to these variables. The 
canonical equations of motion read as 

^-^^ = '2Kzz{q)mz{q) - iKz^{q)ip{q), 

h^^^ = —K^^iqMq) - iKz^{q)mM- (10) 

According to Eq. f llOp the energy of the collective mode with the wave vector q = qi^ 
has the form Vt{q) = K^p^{q)Kzziq) +Kz^{q). Rotating the axes we find the excitation 
spectrum at general q 



r](q) = ^K^^{q)KzM) + K-^M)- (H) 
In what follows we use the reference frame with the x axis directed along Q (Q = (Q, 0)) 
that coincides with the direction of the current jcF Eq- (E]). 

The wave vectors for the stationary waves satisfy the equation 

0(q,i) = 0. (12) 



Stationary waves in a superfluid exciton gas in quantum Hall hilayers 5 



^ 3 














\ 








1 







-0.5 -0.25 0.25 0.5 

V 



Figure 1. Critical interlayer distance versus the filling factor imbalance 

The necessary condition for the observation of stationary waves is the existence of 
nontrivial solutions of Eq. f[T^ . Beside that, the spectrum fl{q) should be real valued 
at all q (complex valued fl{q) signal for the dynamical instability of the state Eq. ([1])). 
These two conditions determine the parameters of the systems where stationary waves 
can emerge. 

The spectrum ( fTTl) depends on three parameters: the ratio d = d/i, the filling factor 
imbalance 9 = {ui — 1^2) /2 and the gradient of the phase Q. The spectrum should be 
real valued at Q = 0. It gives the restriction d < dc{v). The dependence dc{i') is shown 
in Fig. [H The stationary wave equation (JT2ll has nontrivial solutions at Q > Qi{d, z>). 
On the other hand, the condition for the quantity f2(q) be real valued at all q yields 
the restriction Q < Q2{d,u) {Qi < Q2) Thus, stationary waves can be observed at 
Qci < Q < Qc2 that corresponds to jci < 3cf < jc2- The quantities jd and jc2 are the 
function of d and 9. The dependences jd and jc2 on the filling factor imbalance are 
shown in Fig. |2] (these dependences are symmetric with respect to u and we present 
them only for z/ > 0) On can see from Fig. [2] that in balanced bilayers {u = 0) the 
diapason of currents at which stationary waves can be excited shrinks to zero. 

The existence of two critical currents jd and jc2 reflects the general situation 
in superfluid systems. The existence of Landau critical velocity (that in our case 
corresponds to the current jd) is a common feature of superfluid systems. While a single 
component uniform Bose-Einstein condensate does not demonstrate the dynamical 
instability, more complex superflud systems do. The examples are superfluid Bose- 
Einstein condensates in a lattice [281 129], phase-separated two-component superfluid, 
where the components have a common interface [301 El] , and two-component superfluid 
mixtures (with penetrating components) [32 l [331 
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Figure 2. The critical currents jci (solid line) and jc2 (dashed line) at d/l = 0.3 (a), 
d/t = 1.1 (b), and c?/£ — 1.2 (c). The current density j is in e?" jfieP' units. The range 
of currents and filling factor imbalance where stationary waves are excited is shown 
grey. 
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In general case the Landau critical velocities are smaller than ones that cause the 
dynamical instability, but in some cases the system switches directly to the dynamical 
instability regime. As was shown in [32], such a situation takes place in a symmetric two- 
component Bose-Einstein condensate (two components with the same characteristics) 
if the components flow with equal velocities in opposite directions. In the latter case, 
in similarity with a situation in balanced quantum Hall bilayers, there is no range 
of velocities where stationary wave pattern may emerge. Such a similarity can be 
understood as follows. In a quantum Hall bilayer with ut = 1 the exciton transport can 
be interpreted as counter-propagation of two exciton species of different polarization, and 
in balanced bilayers such species move in opposite directions with the same velocities. 

The dynamical instability means that the magnetoexciton state is unstable 
irrespectively to the interaction of magnetoexcitons with an environment. The 
interaction with an environment may result in a Landau (thermodynamic) instability in 
superfluids that develops under the flow with the velocities larger than the Landau 
critical velocities [301 [35] . Phenomenologically, such type of instability can be 
described by a modified Gross-Pitaevskii equation that includes a dissipative term [35] . 
Microscopically, the dissipative term can be connected with an interaction between the 
condensate and normal components in a state where these two component are out of 
equilibrium with each other, and the normal component interacts with an environment 
[361 [37] . As was shown in [311 135], the Landau instability may reveal itself in a vortex 
nucleation at the interface between two superfluid. In a typical situation the Landau 
instability is much weaker than the dynamical instability. In our study we do not take 
into account the Landau instability effects implying that the system is in a regime where 
such effects are small. 

The quantity Kzipiq) in the spectrum flTTl) contains the factor cos^o which sign 
coincides with the sign of the filling factor imbalance. Therefore at z/i > z/2 the stationary 
wave equation ([T2|) can be satisfied for the wave vectors that have negative projection 
on the direction of Q {q^ < 0), and at Ui < z/2, for the wave vectors with q^ > 0. In 
systems that differ only by a sign of the imbalance the stationary wave patterns will be 
the mirror images of each other. 

To be more specific we consider below the case ui > V2- It is convenient to 
parameterize the stationary wave crests by the angle ^ that determines the direction 
of the wave vector counted from the —x axis: q^t = (— gcos^, gsin^). At vi > i>2 the 
Mach cone is situated in the x > half-plane and the Mach angle is om = 7r/2 — ^m, 
where C,m is the maximum angle of deviation of q^t from the —x direction. 

Eq. ( fT2l) can be considered as an implicit definition of the function qstiO- I^ general 
case this equation has several solutions qst,x{C,)i where the index A numerates different 
families of stationary waves. 

To each wave vector q^^ one can put in correspondence the group velocity vector 



The polar angle x that determines the direction of the group velocity Vg relative to the 




(13) 
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X axis is also the function of ^. The dependence x(0 ^^"^ be obtained from Eq. fll3l) 
under substitution q^t = {—qst{.C) cosC,, QstiO sin^- 

A point defect located at the origin emits stationary waves. Their phase at the 
space point r is equal to 



The latter equation gives the dependence r{C,,9). The coordinates of the wave crests 
can be presented in the parametric form 



The parameter ^ belongs to the interval |^| < ^m- The phase takes the values 
6 = ±2ttN + 6in, where = 1,2, . . ., and 6^ is some initial phase that cannot be 
determined in the kinematic approach (for the calculations we put 6in = 0). The sign of 
6 is the same for a given family of waves and can be found from the condition r{C,, 6) > 
that is equivalent to 6'cos[^ + x(0] < 0- 

As was mentioned in Introduction, the location of stationary waves is determined 
by the type of dispersion of the collective mode. In the case considered the dispersion 
depends on the parameter d. Typical dispersion curves at Q = are shown in Fig. [3l 
One can see that at rather small d the dispersion is positive at all q. At larger d the 
dispersion becomes negative at intermediate q. At d close to critical one the spectrum 
contain a minimum at finite q. The wave crests obtained for these three cases are shown 
in Fig. m One can see that in the case of small d the stationary wave pattern is similar 
to one that was predicted and observed in weakly non-ideal superfluid Bose gases, while 
at larger d an additional family of stationary waves emerges inside the Mach cone. The 
minimum in the spectrum reveals itself in an appearance of cusps on the crest lines. 

3. Stationary wave patters in a superfluid magnetoexciton gas 

To compute the amplitudes of stationary waves we will use the dynamical approach [19] . 
Stationary waves emerge as a response on an obstacle that was appeared in the past. 
The casuality principle can be realized by the following choice for the Hamiltonian of 
interaction with the obstacle 



where V^(r) {i = 1,2) is the interaction potential, pj(r) is the electron density operator, 
and r] = +0. We consider the point obstacle V^(r) = t/j5(r) and imply that the 
interaction f|T6l) does not change the sum of local filling factors of the layers. 
In the state ([I]) the energy of the interaction ( fT6i) reads as 



e = qsfr = -qstiOr cos[^ + x(0]- 



(14) 




(15) 




(16) 




I _|_ gjfjQy^^ 



(17) 
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Figure 3. The spectrum of collective excitations (in /el units) at d/l — 0.3 (solid 
line), d/i — 1.1 (dashed line) and d/l — 1.2 (dash-and-dot line) for the filling factor 
imbalance v = 0.25. 



where Uz = Ui — U2- The fluctuating part of f|T7|) expressed in terms of mz{q) has the 
form 

If pit I I piqQy^^ 2^2 

i^/j = ymzjq)^ e-^. (18) 

The interaction f lTSj) modifies the equations of motion: 

r/(^((7) r/ 1 + e^iiQyi'^ 2^2 
= 2ir..(g)m.(g) - ^ir.,(g)^(g) + ^^^^^ ^'^^ 

^ — ci — " ^2^^^^^^'^^^^ ~ iKzv{q)mz{q). (19) 
Taking the partial solution of ( lT9l) at t = 0, we obtain the following expression for 
ruziq) 

m (q) = K,,iq)il + e-^<^y^')e-^ 

AV2^Kzz{q)K^M)-iKzM)~^^vY' 
Using the inverse Fourier transformation for mz{q) 

cos = cos ^0 + 2 y 5^ (g) e*''^ , (21) 

<? 

the expression for the Fourier component of the electron density difference 

Pi(q) -P2(q) = '5,„o$^cos^x^^ e-'^^^-— , (22) 
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Figure 4. Stationary wave crests at i> = 0.25 and (a) — d/£ = 0.3, jcF = 
4.5 • 10-3e3/fe€2. _ ^ ^ -^^ ^ 2.2 • lO-V/fef^. (p) _ ^/^ ^ ^ 2, 

jcF = 1-9 • / heP . The family located outside (inside) the Mach cone is shown 

by dashed (solid) lines. The spatial scale is 1 = ^. 
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and taking the inverse Fourier transformation of 

Ap(r) = pi(r) - p,{v) = \Y1 - P2(q)] e^'^^ (23) 



S 

q 



we obtain 



\n(v\ = V i^^y(g)[l + cos(g^g/^)]e 2 

The first term in Eq. (1241) is the uniform electron density imbalance Apo- The 
nonuniform part p = Ap — Apo describes the density fluctuations caused by the 
stationary waves. Eq. ( 12^ yields the contribution of the modes with the wave vectors 
directed along the x axis. The contribution of modes with general q can be obtained 
by rotation of the coordinate axes. Summing the contribution of all q we find 

^ V ,^q-r -^^^(q)[l + ^"^(lqxQI^')]^"'^ (OK) 

) 2nP 2S ^ ir,,(q)ir^^(q) - (ir,^(q) - ihr^f ' ^ ' 

Replacing the sum with the integral and using the symmetry properties of (!25|) we obtain 

n 1 /"f r°° 
^^'^ = -2^(2^?^^'/^ I (26) 

where 

j( .^ ^ QK^^jq, [1 + cosjqQi'^ sinQJe"^ 

and X is the polar angle for the vector r. 

The integral over q in Eq. fl2Bl) is evaluated from the residue theorem. The 
integration contour is chosen as shown in Fig. [51 The contour 1 in Fig. |5] corresponds 
to the case where the observation point is outside the Mach cone, and the contour 2, 
inside the Mach cone. The integral along the imaginary axis yields the contribution of 
order (r£)~^, while the integral along the real axis - the contribution of order (r£)~^/^. 
We specify the case r£ ^ 1 for which the first contribution can be neglected. In this 
approximation 

where qx{^) are the poles of J(g, ^), and the factor 

fx, cos(e + x)Im(gA(0) <0 . . 

' \ 0, cos(e + x)lm(g,(0)>0 ^ ^ 

indicates whether is the pole q\{C,) inside or outside the contour of integration Fig. [51 
The integral over ^ is evaluated in the stationary phase approximation 

P(') - ^LL^a(6) ^ . (30) 
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Figure 5. The integration contours used for the computation of the integral (|26|) 



Here are the solutions of the stationary phase equation 

^[gA(Ocos(x + 0] = o, (31) 

the function g\{i) is defined as 

^?A(0 = ^bA(Ocos(x + 0], (32) 
and the amphtude Ax{^) is equal to 



^a(0 = -5Aj^^^^^ReSg=,,(5)[J(g,0]- (33) 



One can check that Eq. (l30i) up to the initial phase yields the same wave crests 
as determined by the parametric equations f|T5|) . Indeed, in the limit — )■ the poles 
q\{C) coincide with the stationary wave vectors determined by Eq. ( |T2|) . The stationary 
phase equation (131]) determines the multi- valued function ^i(x) that is reciprocal to the 
function x(^) defined by Eq. (fT3|) . Let us prove the last statement. Considering q and 
^ as independent variables one finds 

oQx oq q 

dil dil ^ dil cos f , , , 

— = —sm^ + — i 34 

oqy oq at, q 
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Taking into account the relation 

dVL 
dq 

and Eqs. ( !T3|) . (!34|) we obtain 



'i=qst{0 



sin ^ — — cos C, 

tanx = . ,%/ . . (35) 

— cos 4 + ^ sm 4 

that yields 

tan(x + = ^ (36) 

Eq. fl36|) is equivalent to Eq. f l3Tl) . Note that Eq. fl36l) gives the explicit expression 
for the function ^{x) that can be used instead of the implicit expression (fT3l) for the 
computation of the wave crest lines. 

The stationary wave density pattern given by Eq. (I30ll is shown in Fig. [6l For the 
computations we use the same parameters as in Fig. Ht. 

4. Stationary density waves in a graphene-dielectric-graphene sandwich 

Stationary waves in quantum Hall bilayers reveal itself in spatially non-uniform local 
electrical fields that emerge due to variations of electron densities in the layers. The 
period of stationary waves is proportional to the magnetic length. For typical magnetic 
fields used for the study of magnetoexciton superfiuidity in GaAs heterostructures 
magnetic length is of order of 10 nm that corresponds to the stationary wave period less 
than 100 nm (see Fig. [6]). For the observation of stationary waves an electrostatic field 
detector should be located close to the electron layer and its spatial resolution should 
be smaller than the stationary wave period. 

In this section we discuss advantages of use bilayer graphene structures for the 
observation of stationary waves in the magnetoexciton gas. The problem of exciton 
superfluidity in double-layer graphene structures was studied theoretically in a number 
of papers [261 EH Ell HQl HH |12]. Here we follow the paper [26], where the theory of 



magnetoexciton superfluidity in graphene bilayers was developed. 

In difference with quantum wells, the graphene layers can be located at the surface 
of the system. It was reported recently on the creation of a bilayer graphene structure 
that has one open layer and separate contacts in each layers ^3]. 

Another advantage of graphene systems is much larger energy distance between 
Landau levels in comparison with one for quantum wells in GaAs heterostructures. 
That allows to use smaller magnetic fields. Smaller field correspond to larger magnetic 
length and larger spatial period of stationary waves in absolute units. 

Let us discuss this question in a more quantitative way. The energies of Landau 
levels in graphene are given by the expression 



£,,(iV)=.g„(iV)5^0^, (37) 
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Figure 6. The stationary wave density pattern for the spectrum with a minimum at 
finite q. The density difference is given in Uz/2tt£'^Ec units, where Ec = jel is the 
Coulomb energy. The spatial scale is 1 = ^. 

where is integer number, and vp = 10^ m/s is the Fermi velocity for graphene. 
Landau levels in graphene have the additional fourfold degeneracy (four components 
correspond to different combinations of spin and valley indices). In a free-standing 
graphene the = level is half-filled, but in difference with the single component 
systems magnetoexciton superfluidity in the graphene bilayers is not possible at half- 
filling [26]. 

The situation is similar to ones for ut = 2 double quantum wells |15] . Electrostatic 
field applied perpendicular to the graphene layers may create an imbalance of filling 
of = Landau level in the layers 1 and 2. We consider the cases ui = 2 + Ua, 
1/2 = 2 — Ua, and z/i = 3 + z/^, z/2 = 1 ~ ^'a (0 < z^a < !)• Here the filling factors are 
defined as the average number of filled quantum states in a given level per the area 27r£^. 
The quantities Uai = and i'a2 = ^ — give the filling factors for the so called active 
component: only for one (active) component the modulus of the order parameter for 
the electron- hole pairing is nonzero: |Aa| = ^y Uail — Va)- 

The eigenfunction for the zero Landau level in graphene coincides with the 
eigenfunction for free electron gas. That is why the dynamics of the active component 
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in the bilayer graphene system is described by basically the same equations that were 
used in the previous section. 

Here we consider the graphene-dielectric-graphene sandwich structure with two 
open graphene layers (two layers separated by a dielectric layer which thickness is much 
larger than the distance between graphenes in graphite). For such a structure the Fourier 
components of the intralayer and interlayer Coulomb interaction read as 

47re2 {e + l)e«'^ + {e - 1)6"'?^ 

Vs[q) - ^ ^^+l)2ggd_(^_l)2g-gd' 

The dielectric constant for the environment Sen = 1 is assumed. 

The empty states in a partially filled Landau level can be considered as holes if the 
Coulomb interaction energy is smaller than the energy distance between Landau levels. 
Let us use the inequality 

Fs{Q) < Ell{1) - Ell{0) = (39) 

as a quantitative condition of smallness of the Coulomb interaction. The quantity Fs{0) 
determines the strength of the intralayer exchange interaction. The right hand part of 
Eq. f l39|l is the distance between the active {N = 0) and the nearest passive (A^ = ±1) 
Landau level for graphene. For the geometry considered in |26j (two graphene layers 
embedded in a dielectric matrix) the quantity -Fs(O) = \/ t: jle^ j et and the inequality 
(!39|) yields the restriction only of the dielectric constant [e > 2). For the graphene- 
dielectric-graphene sandwich the intralayer exchange energy depends on d: 

_^2[(£ + l)e^'^'+(£- l)e--^'] 
^5(0) = — / dxe 2 ^ i i (40) 

and the inequality ( 139|) yields the restriction d > dci{e), where the function dci{e) is 
given implicitly by the equation ^5(0) = ^J2hvpli. On the other hand, the condition 
for the excitation spectrum be real valued yields another restriction on the parameter 
d\ d < dcivaiS), where Va = Va — ^jl is the imbalance for the active component. 
In difference with the situation shown in Fig. [1] the critical distance dc for the 
graphene-dielectric-graphene sandwich structure depends on the dielectric constant. 
The dependences djye) and dc\{e) are shown in Fig. [7l One can see that dc\ < dc 
for e > 4 and under appropriate choice oi d/i one can satisfy the both restriction. At 
fixed d the magnetic length is restricted from above and from below, but the diapason 
of allowed i can be shifted by a change of d. Thus, one can adjust the parameters to 
make the period of stationary waves appropriate for the observation. 

In Fig. [8] we present the stationary wave density pattern computed for a graphene- 
dielectric-graphene structure. The cases of one family and two families without cusps 
are shown. One can see from the Figs. [6] and [8^,b that the amplitude of stationary 
waves grows up when d approaches to dc. 




14 8 12 8 

Figure 7. Critical interlayer distance dc (solid lines 1 - z>a = 0, solid lines 2 
- t'a = 0.25) and dd (dashed line) for the graphene-dielectric-graphene sandwich 
structure. 




0.14 



L-0.14 



Figure 8. Stationary wave patterns in a graphene-dielectric-graphene sandwich 
structure with £ = 15 and = 0.25. (a) - d/i = 1.1, jcF = 1-1 • lO'^e^ /Hei"^; 
(b) - d/i = 1.5, jcF = 6.9 • IQ-^e^/hei'^ 

5. Conclusion 

In conclusion, we wave studied stationary waves in a magnetoexciton gas in quantum 
Hall bilayers. Stationary waves are static excitations of the difference of electron 
densities in the layers. They are induced by the counter-propagating electrical currents 
flowing in a system with an obstacle. It is found that such waves can be excited only 
in imbalanced bilayers in a certain range of currents. The stationary wave pattern 
is modified qualitatively under a variation of the ratio of the interlayer distance to 
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the magnetic length d/i. At small d/£ the pattern is similar to one that emerges 
in a superfluid weakly non-ideal Bose gas: stationary waves are located outside the 
Mach cone. At larger d/£ an additional family of waves is located inside the Mach 
cone appears. Under further increase of d/i cusps are developed in the crests that 
correspond to the second family. The amplitudes of stationary waves grow up under 
increase of d/i. A typical spatial period of stationary waves is several magnetic lengths. 
A convenient system for the observation of stationary waves is a graphene-dielectric- 
graphene structure. In such a system the spatial period of stationary waves can be 
rather large and one can put a detector close to the graphene layer. 

This study is supported by the Ukraine State Program " Nanotechnologies and 
nanomaterials" Project No 1.1.5.21. 
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